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Abstract. We study Bose-Einstein condensation in a linear trap with a dimple 
potential. The present work aims to consider a linear potential in ID decorated 
^ ' with a dimple type potential where we model dimple potentials by Dirac 8 function. 

Attractive and repulsive dimple potentials are taken into account. We also aim to 
compare the results of this paper with the results of a previous paper in which the 
harmonic trap with a dimple potential in ID was investigated [TS]. This model allows 
for simpler, explicit numerical and analytical investigations of noninteracting gases. 
Thus, the Schrodinger equation is used instead of the Gross-Pitaevski equation. We 
calculate the atomic density, the chemical potential, the critical temperature and the 
Q . condensate fraction. The role of the relative depth of the dimple potential with respect 

to the linear trap in large condensate formation at enhanced temperatures is clearly 
revealed. Moreover, we also present a semi-classical method for calculating various 
quantities such as entropy analytically. 
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1. Introduction 

The experimental realization of Bose-Einstein condensates (BEC) of dilute gases of 
alkalis 87 Rb, 23 Na and 7 Li have been observed at very low temperatures following 
the improvement in cooling and trapping techniques [U EJ |3]. Experimentally 
available condensates are systems with finite number of atoms N confined in spatially 
inhomogeneous trapping potentials. 

Several studies jH [5] have revealed that BEC can occur in harmonically trapped 
lower-dimensional systems for finite N despite the enhanced importance of phase 
fluctuations |6J. Bose-Einstein condensation in one dimension (ID) with a harmonic 
trap is attractive due to the enhanced critical temperature and the condensate fraction 
[U Recently, one- and two- dimensional BECs have been created in experiments [7J. 
One-dimensional condensates have also been also generated on a microchip [HI |9] and in 
lithium mixtures [TO] . 

A simple magnetic trap mechanism is called quadrupole trap which varies linearly 
in all directions and vanishes at the center. This trapping mechanism has been used 
to obtain BEC of dilute gases however it has a major disadvantage that there can 
be appreciable losses in the vicinity of the node in the field [11] . There have been 
attempts to overcome this disadvantage [21 [11] of the simple quadrupole trap. Dimple 
type potentials which are used to increase phase space density [121 02] and condensate 
fraction may be useful to circumvent or the loss problem of quadrupole trap. The 
Schrodinger equation equation for a particle in the V-shaped potential decorated by a 
repulsive or attractive Dirac delta function interaction at the center has been solved very 
recently [H]. The present work aims to consider a linear potential in ID decorated with a 
dimple type potential where we model dimple potentials by Dirac 5 function. Attractive 
and repulsive dimple potentials are taken into account. We also aim to compare the 
results of this paper with the results of a previous paper in which the harmonic trap 
with a dimple potential in ID was investigated [15]. On one hand, the linear trap 
could be imagined as a disadvantage against harmonic trap applications. However, our 
investigations do aid to comprehend the conditions of efficient loading of atoms in the 
trapping techniques. We calculate the transition temperature as well as the chemical 
potential and condensate fraction for various number of atoms and for various relative 
depths of the dimple potential. For describing a system with interacting particles, the 
Gross-Pitaevski equation is usually utilized. We note that we neglect the interactions 
between the atoms in our model, and thus the Schrodinger equation for the linear trap 
with the dimple potential is solved. 

Cavalcanti et.al. showed even 3D gases in the presence of a uniform field do not go 
undergo Bose-Einstein condensation at finite temperatures [16]. However, the authors 
have revealed that if there is a point-like impurity at the bottom of the vessel, Bose- 
Einstein condensation occurs at T ^ 0. Moreover, as mentioned in a previous paper 
[15] . modification of the shape of the trapping potential can be used to increase the 
phase space density [12]. "Dimple" -type potentials are the most favorable potentials for 
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this purpose. The phase-space density can be enhanced by an arbitrary factor by using 
a small dimple potential at the equilibrium point of the harmonic trapping potential 
[T7] . The effect of variations of the longitudinal trap frequency on the BEC profile 
has been analytically treated by using delta kicks [18]. Ground state properties of ID 
quantum gases in a harmonic trap in the presence of a point-like attractive potential 
have been studied very recently [19]. The formation of a Bose-Einstein condensate in a 
cigar-shaped three-dimensional harmonic trap, induced by the controlled addition of an 
attractive "dimple" potential along the z axis has been studied very recently [20]. We 
model the dimple potential using a Dirac 5 function. The delta function can be defined 
via a Gaussian function [21] g(x,a) = (l/y^a) exp (—x 2 /a 2 ) of infinitely narrow width 
a so that g(x, a) — > 5(x) for a — > 0. This allows for analytical calculations in some 
limiting cases as well as a simpler numerical treatment for arbitrary parameters. 

The paper is organized as follows. In Sect. II, we present the analytical solutions of 
the Schrodinger equation for a Dirac 5-decorated linear potential and the corresponding 
eigenvalue equation. In Sect. Ill, we give semiclassical approach and calculate the 
entropy of the system. In Sec. IV, determining the eigenvalues numerically, we show 
the effect of the dimple potential on the condensate fraction, the transition temperature 
and the chemical potential. In the limit of strong dimple potential the results are 
presented in Sect. IV. Finally, we present our conclusions in Sect. V. 

2. Linear confining potential decorated with a Dirac 5 potential at the 



We begin our discussion with the one dimensional linear potential decorated with a 
Dirac 5 function at x=0, {S(x)) . This potential is given as: 



where m is the mass of the particles, / is force term due to the linear trap and a is the 
strength (depth) of the dimple potentials located at x — 0. The factor ft 2 /2m is used for 
calculational convenience. The negative a values represent repulsive interaction while 
the positive a values represent attractive interaction. The time-independent Schrodinger 
equation for the potential given in Equation (OQ) is: 



Although the solution of this equation can be found easily we summarize it for the 
sake of completeness. First, we solve Equation (EJ) for x =fi 0. Since the potential is 
an even function the eigenstates of this potential are even and odd functions (see e.g. 
|22j). Hence it is enough to find the solution for positive x. Defining k 2 = 2mE/h 2 , 
I = (h 2 /(2mf)) 1/3 and u = x — k 2 l 2 and inserting these definitions into Equation ([2]) 
one gets the Airy Differential equation for the defined variable [23] : 
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Therefore for o = in Equation (j2J), utilizing the fact that the wave functions are 
odd and even one finds the wave functions in terms of Airy Ai function Ai(x) as [23] 
(Since the second solution of AiryBi function diverges for x — > oo the wave functions 
are expressed only in terms of Airy Ai function): 

^ e (x) = < Ai[~(\x\ - k 2 J 3 )} for n = 0, 2, . . . 

* (x) = c n sgn(x)Ai[j(\x\ - k 2 n f)} for n = 1, 3, . . . . (4) 

Here, c n and c° are the normalization constants for even and odd wave functions 
respectively and sgn(x) stands for the sign function. The c e n and c° can be found 
using the roots of Ai(x) and the derivative of AiryAi function Ai'(x) [23]. Both Ai(x) 
and Ai'(x) have infinite number of zeros and all the roots of these functions are negative. 
Let us denote the roots of Ai(x) by — s n and the roots of Ai'(x) by — z n , i.e Ai(—s n ), 
Ai'(—z n ) = 0. These coefficients can be ordered as —z\ < —s\ < —z 2 < — s 2 < . . . 
. Ordering these coefficients just using one symbol like a = Z\,a\ = Si,...a 2n -i = 
s n , «2n-2 = z n , n = 1,2, . . ., one gets for the normalization coefficients [23] : 

c e n = [2a n l[Ai{-a n )fY l/2 for n = 0,2,... 

c° n = [2Mi'(-a n )] _1 for n = 1, 3, . . . . (5) 

Using the fact that ^ o (0) = and (cNf e (x)/dx) x= o = 0, one finds the eigenvalues of the 
Hamiltonian for the potential given in Equation fll]): 

We define 



E n = ( -t— ) a n for n = 0, 1, . . . . (6) 



k n = y/(2mE n )/h. (7) 

Using I = (h 2 /(2mf)) 1 ^ 3 we get 

k n = yfotn~/l for n = 0, 1, . . . . (8) 

If cr in Equation (j2J) is not zero the eigenvalues of the even states can be found by 
the taking the integral of the Schrodinger equation from — e to e in the limit e — > as 
outlined in [21] [25]: 

2A»'(-Aff 2 ) 

which can be solved numerically for k n . Substituting these k n values into the argument of 
AiryAi function for even wave functions in Equation (j3J) one gets the even eigenfunctions 
for cr 7^ of the Hamiltonian for the potential given in Equation (pQ). The odd 
eigenfunctions and their eigenvalues are unaffected because the Dirac 5 function is 
located at x = and the odd eigenfunctions vanish at this point. On the other hand, 
we showed that the energy eigenvalues of even states change as a function of a. The 
ground state energy eigenvalue decreases unlimitedly as a increases (attractive case). 
However, the energies of the excited even states are limited by the energies E 2n +\ of 
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odd states and as a — > oo, E 2n+ 2 — > E 2n+ i = (2n + 1 + 1/2) fio; where n = 0,1,.... 
Thus, as a — > oo the energy eigenvalues of these states go to the energy eigenvalue of 
the one lower eigenstate, and the odd energy eigenstates asymptotically become doubly 
degenerate. 



3. Semiclassical Approach 

We will first approach the problem semiclassically. The spectrum of the Hamiltonian for 
the potential given in Equation ([1]) is discrete for all values of a including o = (only 
the linear confining potential). For convenience, we define a dimensionless parameter A 
in terms of a as: 

A = al. (10) 

We choose the value of I such that we get the same order for A values for a given a 
as in Ref. [15], so that one can compare the effect of dimple potential to the linear 
confining potential with the effect of dimple potential to the the harmonic trap. In [15] 
the calculations are made for 23 Na particles with m = 23 amu and the frequency u of the 
harmonic trap was chosen as u = 2n2lHz. Equating the length scales of harmonic trap 
yjh/ (2muj) to the length scale of the confining linear potential /, we get for I = 3.23/xm. 

Quantum mechanically - as we will present in the following section- one calculates 
the thermodynamic variables of a non-interacting condensate like critical temperature, 
condensate fraction and the chemical potential taking sums over the discrete energy 
eigenvalues. However, in the semiclassical approximation one uses the density of states 
and employs integrals instead of sums (see e.g. [TT1 [26]). Therefore we first find the 
density of states for a linear confining potential. As we have showed in the previous 
section, the odd and even energy eigenvalues of the linear confining potential (a = in 
Equation (CQ)) is proportional to the negative of the root values of Ai(x) and to the Ai'(x) 
(s n and z n values defined in the previous section), respectively. This proportionality can 
be used to calculate the density of states for linear confinig potential. The s n values can 
be approximated by the formula [23] , [j]: 

sn = [^{n-\)fl\ (11) 

We find the density of states using s n = {2m/h 2 f 2 ) l ^E n and assuming n and E n 
as continuous variables. Doing this we get 

dE~ A 2 / 2J ' 1 ] 

However this is only the density of odd states. For each odd state there is an even state. 
Hence the density of states (g(E)) for the linear confining potential is 

= I < JP>)" 2 £ I/2 - (") 

| This approximation can be improved to get a better accuracy |23) . However, since we are only 
interested in the density of states and more accurate formulae give the same density of states we use 
the simplest approximation to the roots. 
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Using the density of states one can calculate the number of particles in the 
thermal gas that is the number of particles which are not in condensate phase [llj. Above 
the critical temperature T c the number of particles in the ground state is negligible and 
all the particles are assumed to be in the thermal gas. As the temperature decreases the 
chemical potential increases and comes very close to the ground state energy. Therefore 
one gets the critical temperature taking the chemical potential \i is equal to the ground 
state value [26jjj|. Using the density of states and Bose distribution we find: 



ei 3 c (e-A0-l ^fitjl' e p c (e-n)-l 

where the lower limit of the integral is the energy eigenvalue of the first excited state 
and /3 C = 1/(/sbT c ). For different A values which are used to model the strength of the 
dimple potential we can find the energy eigenvalue of the ground state using Equation 
((9]). The energy eigenvalue of the first excited state is same for all A and given by 
Equation (J6]) taking n = 1 because it is an odd eigenf unction. The Dirac 5 potential in 
Equation (CQ) does not change the density of states. That is it does not create or destroy 
eigenstates but only change the values of the energy eigenvalues. Therefore the same 
density of state expression can be used for all A values. However, the total number of 
particles in the thermal gas or (in the condensate) differ for different A values because 
as A varies the energy value of the ground state varies which cause a change in the 
chemical potential value. 

So, if the Equation (|14p can be solved for /3 C for fixed N, one can find the change 
of critical temperature with increasing o which determines the strength of the dimple 
potential. We approximately calculate the integral in Equation (I14p taking the lower 
limit zero (i.e we set ex = 0) and shifting the chemical potential by the same amount 
(n — > /i — ei). Applying these approximations and taking the integral in Equation (fl4l) 

^-2(-|^) 1 / 2 (i-)3/ 2 ,3 /2 ( e ^). (is) 
where the Bose function g p (x) is defined as [TT] : 

oo i 

9*V) = E r P ■ (16) 

1=1 

We have solved this equation and find T c for various A values. We will present our results 
in the following section in order to compare the results obtained from semiclassical 
approximation with the results obtained from quantum mechanical calculations. 

Adding a dimple potential adiabatically changes the critical temperature by 
changing the value of chemical potential fi [26J. If an attractive dimple is added to 
the trap potential at temperatures slightly above the critical temperature of the trap 
without a dimple adiabatically, the entropy remains constant. Hence the decrease in the 
ground state energy cause the formation of BEC [26]. In the laboratory this process can 

§ As mentioned in Ref. [4] , the phase transitions due to discontinuity in an observable macro parameter 
occurs only in thermodynamic limit, where N — ¥ oo. Therefore we assume the critical temperature T c 
is the temperature that particles begin to accumulate in the ground state. 
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be maintained adiabatically and the entropy of the boson gas remains constant during 
the process. Therefore if one finds an expression for the entropy of the Bose gas one can 
calculate the value of condensate fraction of a BEC obtained for an adiabatic process. 

The entropy of a non-interacting Bose gas is easily calculated using the grand 
potential Q [26J: 

n = Q + k B Tj2 ln l 1 -e' 3{fM ~ et) }- (17) 

i 

Using the density of states for a linear confining potential given in Equation (fl3|) . the 
grand potential of a non-interacting Bose gas trapped using linear confining potential 
can be written as 



n = n + ^(^) 1 / 2 r ln[l - e^] de. 



(18) 



Calculating the integral in Equation ( TT8|) one finds 

2 2m 



/3 2 J + 2/35/2 ^ n 5 /2 



(19) 



where (3 = l/(k B T), Erfc(x) = 1 — Erf(x) = 1 — J Q x e t2 dt and g p (x) is defined in ffT6]) . 
Then, we find the entropy by S = —dQ/dT: 

{^92(e^-^) + yfcQi - ei ) ln[l - e^ 1 )] + 

5^F e w ^Erfc(Vh^e7) y^x y> e n/ ^Erfc(y/h^) ) 

n=l n=l 

(20) 

For the temperatures slightly above T c , one can take /x = eo and /3 = /5 C . Therefore for 
temperatures slightly above T c , Equation( l20~i) becomes 

{^^ 2(e/3c(eo ~ ei)) + ^ (e ° " ei) ln[1 " e/3c(£0 " ei)1 



2/3, 



^ / e n/3ce °Erfc(Vn^e7; 



U V" \ ^ 



o /p3/2 Z— / 7^5/2 
z Pc n=1 



y^p e^ £ °Erfc(Vh^7) 
2/3 c l/2 £j n3/2 



(21) 

where eo = (^~) 1//3 ) #i an d ei = (i^ - ) 1 ^ 3 ) s i- Adding a dimple potential changes the 
ground state therefore the chemical potential changes also due to \i = e g < e . So, we 
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Figure 1. The condensate fraction values evaluated with the semiclassical method for 
N = 10 4 . The logarithmic scale is used for x axis, where AE = €q — e g . 
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S = k R - 



2 , 2m 



,1/2 
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^( es _ ei )l n [l_ e A(M- £1 )] 



3/2 



2/3 c l/2 



D 

n=l 

oo 



e n ^ £ »Erfc(v^cii) 



n 5/2 



n 



3/2 



(22) 



for the entropy of the BEC in the linear trap with a dimple potential. Since we model 
the dimple potential with a Dirac 5 potential, we can calculate the new ground state 
energy (e s ), -i.e the ground state energy of the trap with a dimple- using Equation fl9]). 



Then, we equate Equation fl2Tl) and Equation f[2"2l to calculate the final temperature of 
the system, because the entropy remains constant during adiabatic process. Then, we 
calculate condensate fraction for this final temperature. 

We present the increase of the condensate fraction with the decrease in the ground 
state energy value in Fig. (OQ). To obtain this figure, we have first calculated the 
ground state energy values (e 9 's) of linear confining potential with a dimple potential 
for different strengths of the dimple i.e for different a (A) values in our model. We choose 
the range of a values which determines the strength of the dimple same as in the paper 
|15j . i.e. between Om" 1 and 10 8 m _1 . Finally we used Equation f|T5|) to find the number 
of particles in the thermal gas Nt- Since we assume that the number of particles N in 
the Bose gas is fixed we could get the condensate fraction using (N /N = 1 — N T /N). 
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Figure 2. (a) The critical temperature T c vs. A for N — 10 4 and N = 10 5 and for 
negative A. A is a dimensionless variable defined in Equation (fTOj) . The solid line 
and the dashed line show T c values with decraeasing A for N = 10 4 and N = 10 5 , 
respectively. Here we use m = 23 amu ( 23 Na) and I = 3.23/im. The logarithmic scale 
is used for —A axis, (b) The critical temperature T c vs. A for N = 10 4 and N = 10 5 
and for positiveA. The logarithmic scale is used for A axis. 



4. BEC In a One-Dimensional Linear Potential with a Dirac 5 Function 

In this section we will investigate the behavior of a BEC confined in a one-dimensional 
linear confining potential decorated with a Dirac 5 function given in Equation ([2]) 
quantum mechanically. That is we will use the discrete energy eigenvalues calculated 
numerically by Equation ([8} and Equation fl9]) for odd (n = 1, 3, • • •) and even 
(n = 0, 2, ■ • •) states, respectively. 

First we will investigate the change of the critical temperature as a function of A 
which shows the effect of the dimple potential in our model. We will present the change 
of the critical temperature with respect to A both for attractive and repulsive dimple. 
The critical temperature (T c ) is obtained by taking the chemical potential equal to the 
ground state energy (// = E g — E ) and 

i=l 

at T = T c , where j3 c = 1/(&bT c ). For finite N values, we define T c ° as the solution of 
Equation (123 j) for A = i.e. only for the linear confining potential. 

In Equation (123]) . £j's are the eigenvalues for the potential given in Equation (JT|). 
For A / the energies of odd states are unchanged and equal to s/a^/l for n — 1,2, 3. 
The energies of even states are found by solving Equation Q numerically. Then, these 
values are substituted into Equation (123]) ; and finally this equation is solved numerically 
to find T c . 

For positive A, we present the change of T c with increasing A for = 10 4 and 
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Figure 3. The critical temperature T c vs. A for N — 10 4 for harmonic and linear trap 
with a dimple, respectively. A is defined in Equation (fpQ|). The solid line and the dashed 
line show T c values with increasing A for harmonic and linear traps, respectively. The 
other parameters are same as Fig. ([2]) The logarithmic scale is used for A axis. 

N = 10 5 in Fig. ([2]-b). For negative A, we present the change with increasing absolute 
A in Figure pj-a). In this figure, logarithmic scale is used for |A| axes. As |A| increases, 
the critical temperature increases (decreases) very rapidly when |A| ~ 10 2 for positive 
(negative) |A|. However the decrease of critical temperature for negative A is not as 
rapid as the increase for positive A. This is due to the fact that as |A| increases for the 
attractive case the ground state energy value decreases indefinitely however the excited 
states are bounded. The limit value of the eigenenergies of the excited states for with 
a tough dimple (modeled by a Dirac 5) are the eigenvalues of the one lower level of the 
confining potential without a dimple. On the other hand for the repulsive case (negative 
A in our case) all the energies with a tough dimple are bounded above by the eigenvalues 
of one higher level of the confining potential without a dimple. 

At this point it may be useful to compare the critical temperature values of a 
harmonic trap and linear trap with in ID with the increasing depth of dimple (A values 
in our model) [jj]. We show the change of the critical temperature with respect to A for 
both cases in Fig. ([3]). The critical temperature values are higher for harmonic trap 
compared to linear trap if one chooses the length scales same for both potential. This 
is an expected result because the confinement strength of harmonic potential which 
increases quadratically with respect to distance to the center of the potential is larger 
than linear trap which increases linearly with respect to distance to the center. 

As we mentioned in the previous section the increase of the critical temperature 
can be calculated approximately by solving Equation ( TH|) obtained by the semi-classical 
method. We compare the critical temperature values obtained quantum mechanically 

|| The change of the critical temperature with respect to increasing depth of dimple for a harmonic 
trap is presented in Ref. [15] . 
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Figure 4. (T c sc ', T c ) vs. A for N = 10 4 . The dashed line shows the T c sc values 
obtained semiclassically where the solid line shows the T c values obtained quantum 
mechanically. The logarithmic scale is used for A axis. 



and semiclassically in Figure (Hj). This figure shows that the agreement of the results 
of two methods is better for large A (or large AE) values. This is due to the fact that 
the average occupations for excited states is so low (n excited <^ l) that the semiclassical 
treatment for this gas is valid. Therefore, the semiclassical method may be useful for 
calculations with large A values because for A »s 10 3 the numerical calculations do not 
give accurate results. 

For a gas of N identical bosons, the chemical potential /x is obtained by solving 
equation 

N = V S i = N + V -37 r , (24) 

i=0 i=l 

at constant temperature and for given N, where Si is the energy of state %. We present 
the change of /x as a function of T/T c ° for iV = 10 4 in Fig.flH]) and in Fig. ([6]) for A = 
and A = 32, respectively. By inserting /x values into the equation 

N « = e ^-!.)-l - < 25 ' 
we find the average number of particle in the ground state. N /N versus T/T c ° for 
= 10 4 , A = 0, 3.2, 32 are shown in Fig. (J7|). We also present the change of condensate 
fraction for iV = 10 4 and iV = 10 5 with increasing temperature for fixed A = 32 in 
Fig. (151) . The condensate fractions for different N values are drawn by using their 
corresponding T c ° values. The T c ° values are 0.37/xK for iV = 10 4 , and 1.7/xK for 
N = 10 5 . 

We also find the condensate fraction as a function of A at a constant T = T c °. These 
results for iV = 10 4 are shown in Fig. ([9]). We notice that large A values (A > 1) induce 
sharp increase in condensate fraction like in the harmonic trap with a dimple potential 

im 
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Figure 5. The chemical potential /i vs temperature T/T® for N — 10 4 and A = 0. 
The other parameters are the same as Fig. [5] 




0.5 1 1.5 2 2.5 3 



TH° C 

Figure 6. The chemical potential /i vs temperature T/T® for N = 10 4 and A = 32. 
The other parameters are the same as Fig. [5] 

Finally, we compare density profiles of condensates for a linear trap and a linear 
trap decorated with a delta function (A = 4.6) in Fig. (fTUjl . Since the ground state 
wave functions can be calculated analytically for both cases, we find the density profiles 
by taking the absolute square of the ground state wave functions. 

5. Conclusion 

We have investigated the effect of the tight dimple potential on the linearly confined 
one dimensional BEC. We model the dimple potential with the Dirac 5 function. 
In our model system, the increase in the strength of the Dirac 5 function can be 
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0.5 1 1.5 2 2.5 3 



T/T o ° 

Figure 7. N /N vs T/T° for N = 10 4 dotted line for A = 0, dashed line for A = 3.2, 
solid line for A = 32. The other parameters are the same as Fig. [2j 




0.5 1 1.5 2 2.5 3 



Figure 8. Condensate fraction N Q /N vs T/T c ° for fixed A = 32. The solid and dashed 
lines show condensate fraction for N = 10 4 and N = 10 5 , respectively. Logarithmic 
scale is used for A axis. The other parameters are the same as Fig. O 

interpreted as increasing the depth of a dimple potential. This allows for analytical 
expressions for the eigenf unctions of the system and a simple eigenvalue equations 
greatly simplifying numerical treatment. We have calculated the critical temperature, 
the chemical potential and the condensate fraction and demonstrated the effect of the 
dimple potential. We have found that the critical temperature can be enhanced by 
an order of magnitude for experimentally accessible dimple potential parameters. In 
general, we find that T c increases with the relative strength for a non-interacting gas 
for an attractive dimple and decreases for repulsive dimple. However as our results 
show the decrease of the critical temperature for the negative dimple is not as rapid 
as the increase for the repulsive case. These results may be used to propose that the 
dimple potentials can be a good candidate to circumvent the disadvantage of the simple 
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A 

Figure 9. Condensate fraction No/N vs the strength of the Dirac S potential A for 
N = 10 at T = T®. The logarithmic scale is used for A axis. The other parameters 
are the same as Fig. [2j 




Figure 10. Comparison of density profiles of a BEC in a harmonic trap with a BEC 
in a harmonic trap decorated with a delta function (A = 4.6). The solid curve is the 
density profile of the BEC in decorated potential. The dashed curve is the density 
profile of the ID harmonic trap (A = 0). The parameter z is dimensionless length 
defined after Equation ([2]). The other parameters are the same as Fig. [2j 



quadrupole trap mentioned in Reference jTT]. Repulsive dimple can be used to repel 
atoms from the vicinity of the node (x = 0) of the quadrupole trap. Moreover, It may be 
also possible to use attractive dimple potentials for the same reason because only dimple 
type potentials are also able to trap several atoms [13J. So they can prevent loosing 
atoms in the vicinity of the node. Moreover the increase of the critical temperature will 
be a further advantage for attractive dimple potentials. 

We have also analyzed the change of the condensate fraction with respect to the 
strength of the Dirac 5 function at a constant temperature (T = T c °), and with respect to 
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temperature at a constant strength. It has been shown that the condensate fraction can 
be increased considerably and large condensates can be achieved at higher temperatures 
due to the strong localization effect of the dimple potential. Finally, we have determined 
and compared the density profiles of the linear trap and the decorated trap with the 
Dirac 5 function at the equilibrium point using analytical solutions of the model system. 
Comparing the graphics of density profiles, we see that a dimple potential maintain a 
considerably higher density at the center of the linear trap. 

The comparison of the critical temperature values for BEC s in linear trap with a 
dimple and harmonic potential with a dimple shows that the critical temperature values 
are higher for harmonic trap for the same dimple strength which is an expected result 
since the confinement of harmonic potential are more powerful. However the simplicity 
of the linear trap may still be useful in obtaining BEC. 

We have also presented a semi-classical method for calculating various quantities 
such as entropy, critical temperature etc. The results show that as the dimple strength 
increases the semi-classical approximation gives better results for non-interacting gas. 

We believe that the presented results which are obtained for the case of 
noninteracting condensate by using delta function modeling of tight dimple potentials 
can help significantly such experimental examinations. 
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